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Abstract - In the present paper, we consider a model of non-minimal modified Yang-Mills 
(Y-M) theory in the Friedmann-Robertson- Walker (FRW) cosmology, in which the Y-M field 
couples to the scalar curvature through a function of its first invariant. We show that cosmic 
acceleration can be realized due to non-minimal gravitational coupling of the modified Y-M 
theory. Besides general study, we consider in detail the case of power-law coupling function. 
We derive the basic equations for the cosmic scale factor in our model, and provide several 
examples of their solutions. 
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1 Introduction 

Present accelerated expansion of the universe is well proved in many papers [l]-[8]. In order 
to explain so unexpected behavior of our universe, one can modify the gravitational theory [5]- 
[14], or construct various field models of so-called dark energy which equation of state satisfies 
^SJ ■ 7 = p/p < —1/3. The most studied models consider a canonical scalar field (quintessence) [15]- 

[T7] . a phantom field, that is a scalar field with a negative sign of the kinetic term [TH]-[1I], or the 
OO . combination of quintessence and phantom in a unified model named quintom [22)-|25|. In such 

field dark energy scenarios, the potential choice plays a central role in the determination of the 
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lO . cosmological evolution. 

' The alternative approach to the problem of accelerated expansion is the consideration of various 

modifications of the gravity theory. Among such modifications, the theories with non-minimal 
coupling of a field to gravity are especially attractive [26]- [32]. For instance, scalar tensor theories 
are generalization of the minimally coupled scalar field theories in a sense that here the scalar 
I field is non-minimally coupled with the gravity sector of the action i.e with the Ricci scalar R. 

^ • In non-minimal theories, the matter field participates in the gravitational interaction, unlike its 

counterpart in the minimally coupled case where it behaves as a non gravitational source ^33]- [36). 

In particular, the discovery of the accelerated expansion of the Universe encourages further 
development of the Y-M theory, including its non-minimal coupling to gravity. Numerous attempts 
have been made to consider modified Y-M theories in cosmology as alternatives to the dark energy 
(see, e.g., [37] and references therein). One of the directions in such a generalization of the Y- 
M theory is connected with a non-minimal extension of the Y-M field theory. There are several 
motivations to study non-minimal Y-M theory (see, e.g., [3H||31] and references therein). 

Basically, the present study is a sequel to the paper [JU], where the modified Y-M theory 
is investigated under the condition of minimal coupling to gravity. Following the considerations 
provided in [3D], we study a model of non-minimal modified Y-M theory, in which the Y-M field 
couples to a function of the scalar curvature. We show that the accelerated expansion can be 
realized due to the non-minimal gravitational coupling of Y-M field. Besides general study, we give 
some examples of exact solution for the model under consideration, which surely can not cover all 
possible applications of this research. 
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2 The model equations 

Let the action of our model be presented by a generalization of the modified Y-M action [40] on 
the case of non-minimal coupling: 



(1) 



where the Yang-Mills tensor is i^^;^ = diWI^ - dkW°- + f°''"'W°W^, and * and $ are the arbitrary 
differentiable functions. In the case of 4* = 1, $ = F-iF""^^ , this action describes the Einstein- 

Yang-Mills theory. We assume that the universe space-time is described by a Friedmann- Robertson- 
Walker (FRW) geometry: 



(2) 



where ^(r) = sinr, r, sinhr in accordance with a sign of the curvature k — -1-1,0,— 1. To study the 
model based on action ((l} in metrics ©, we substitute this metrics into action ([T]) and take into 
account that 

a'dN - adN + a^N 



R = -6- 



,2^3 



for this metrics. As a result, we can obtain the following effective Lagrangian per unit solid angle 

3 /a^a aa^ a?dN 



/a a aa 
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(3) 



where / = p^^p"-'-^ _ At the same time, the generalized Wu-Yang ansatz for 5*03 Yang-Mills field 
can be written down as 1411: 



= X' 



Substituting 



K{r,t) = P{r) cos a{t), S{r,t) = P{r) sin a{t), W{r,t)^a{t) 
into this ansatz, we get the following components of Y-M tensor [41] 
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Foi Fo2 Fo3 = 0, Fi2 — e P'(r)(^m cosa + 1 sina 
'^P'{r) sinfl^m sina - 1 cosa^ , F23 = e"^ sm6(^P^{r) - l^n 



(4) 



presented in the orthonormalized isoframe n = (sin cos 0, sin0 sin cos 0), 1 = 
(cos cos (/), cos sin — sin 0) and m = (— sin0, cos </>, 0). The prime in Q stands for the deriva- 
tive with respect to r. As it has been noted in ^T], Y-M field (jH) possesses only magnetic compo- 
nents. From formulas and it is easy to find that the Y-M invariant / = pa^paik ^^^^ ^^iq 
following expression: 



2P' 
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Varying the Lagrangian density ([3]) over P(r) and taking into account (O, we obtain the following 
Euler-Lagrange equation for the Y-M field: 



-1 r n 3T T 

\ Q{tW + a^<i>' +P'— Q(m" + a^<i> 
J I i or I 



0, 



(6) 



where = d^{l)/dl, = d'^{l)/dl and the following notation is temporarily introduced: 



Q{t) 



3 /a-'d ad a dN 



+ kaN^ . 



The particular solution of the Y-M equation in the FRW metrics has been obtained in [JT] . It 
has the form P{r) = ^'(r) = cos r, cosh r for k — +1,-1 consequently. It satisfies equation (|6]) as 
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well. Indeed, it turns the expression in braces into zero and satisfies, in view of ([S]), the equality 
dl . 

— = 0. It is easy to prove the latter, as for this solution the Y-M invariant depends only on time; 
or 



I = I{t) 
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e2a4(t)- 



(7) 



Varying Lagrangian over a{t) and N{t) with the subsequent choice of gauge iV = 1, one can obtain 
the following equations for our model: 



2-+ - +— */ -4/*'/ 2--2 - +— +16 - /^^r"/ 

a \aJ a^J La \aJ a^J \aJ 
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Several examples of exact solution 
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(I) The arbitrariness of differentiable functions and $(/) essentially complicates the general 
analysis of the model equations ([S]),®- Therefore, we will consider some special cases. Let us begin 
with the power-law dependencies: 



$(/)=A/", ^{I) = Br 



(10) 



where A and B are some dimensional constants, free parameters of the model. The substitution 
of expressions (flUl) into equations © leads to the following set of equations: 



(1 - 4m) 



2^ 
. a 



- (l-4m) + - 



- (l-4m) + - 



^|(3-4n)r 
3 B 



(11) 
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As one can see from these equations, the case m = 1/4 is a specific one. Indeed, from equation ([TT|) 
with m — 1/4 it immediately follows that n — 3/4, i.e. the power of / in equation is equal 
to n — m = 1/2 , and the first term in the left-hand-side of this equation is equal to zero. Taking 
into account (O, we can reduce equation ([T^ to a ratio between constants A and B of the form: 

8V67rG ^ , , , . . , . , . 

k = . ihus, the scale factor remams uncertam, or arbitrary. It is easy to prove that 

3e B 

the same unclearly interpreted result will turn out, if functions (|10p are directly substituted into 
Lagrangian Q together with n ^ 3/4, m = 1/4 (or $ = A{6/e'^f/'^a-^{t), = A{6/e'^)^^'^a''\t)). 
Varying then it over a{t) and N{t), we again arrive at the same relationship between A and B. As 
for the equation for the second derivative of the scale factor, it will be satisfied identically for any 
a{t). Putting this case aside, we consider our model with m 7^ 1/4. 
Combining ([TT|) and ([T^ . we have the following equation: 



d ^ SttG A [2{n - m) - 1] ^„_, 
a 3 B (4to-1) 



(13) 



for the second derivative of the scale factor. 

It follows from ([7]) that J > 0. Therefore, the necessary condition of accelerated expansion of 
the universe (a > 0), following from equation ()13p . is reduced to 



A [2(r 



.)~1] 



B (4to-1) 



> 0. 



The latter inequality must be solved for m > 1/4, or m < 1/4. The result of solving is presented 
in Table 1. 
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la 


A/B > 


ui ^ 1/4 


n - — 777 \ 1 /O 

/ t; lit ^ A- 1 Lj 


lb 


A/B > 


TO < 1/4 


n — m < 1/2 


Ic 


A/B < 


TO > 1/4 


71 - m < 1/2 


Id 


A/B < 


TO < 1/4 


n — m > 1/2 



Table 1. Conditions for the free parameters A/B, m and n — m, corresponding to 
corresponding to accelerated expansion according to equation ((T3)) . 



It is necessary to emphasize that equation ([T3| is employed for the analysis of the accelerated 
mode of evolution, but is a differential consequence of equation ([T2l) . Therefore, the model dynamics 
is defined by the only independent equation ([T2l) . As one aims to solve this equation at the condition 
of accelerated expansion, it is necessary to take into account Table 1. 

From Table 1, it follows that in the case of standard Y-M Lagrangian {n — 1) the accelerated 
expansion is possible, if (la) A/B > and to £ (1/4,1/2), or (Id) A/B < and to < 1/4 
are realized. The assumption of minimal coupling (that is m = 0) corresponds to (Id). However, 
equation p2|) at A/B < 0,m = 0,n — 1 and (fc = +1) has no the real solution, and the case (la) 
does not correspond to to = 0. The latter simply means the absence of the accelerated mode in 
Einstein- Yang- Mills minimal theory [37]. The trivial solution with accelerated expansion can be 
obtained in the case A/B < 0, to = 0,n = 1, and for the negative sign of curvature. Taking into 
account ([7]), and supposing A — l/167r, B — —1, we can express equation ([T2|) and its solution as 
follows: 

G 1 _ [G 



SttG 
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Ve2 j 


B 



a -1--— , a(t) = (14) 
e a V e 

where the constant of integration is equal to zero for the sake of simplicity. For this solution, the 
acceleration equals d = (G/e^)/[(G/e^) + t^f^^^ > 0, and it decreases to zero with time. It is 
interesting that this model of non-massive Y-M field with linear dependance on / = pa^pnk jgg^^^g 
to the accelerated expansion. However, due to equation (jl2p . this model violates the weak energy 
condition. Therefore this example of solution (|14p is exclusively illustrative. 

Let us consider now two examples of exact solution for this model which are not so trivial but 
rather simple. 

a) Let n = l, TO==3/4, that is we will consider the case (Ic) from Table 1. Then equation ([13 
can be written down as follows: 

2 — k ~ C a, where G — 

The obvious solution for this equation is: 

a(i) = |(t_Go)^-^, 

where Gq > ^k/C"^. This model experiences constant acceleration: d — G/4, and the Hubble 

b) We now consider an example of exact solution with n — m that corresponds to (lb) and (Ic) 
in Table 1. In these cases, equation (IT^ can be written down as follows: 

Md^-5k^Ca^, (15) 

where 

M=|4to-1|, 6 = 



G = ^ 



+ 1 for TO > 1/4 , 
-1 for TO < 1/4 . 



For the positive sing of curvature (fc = -1-1) and to > 1/4, as well as for the negative sign of 
curvature {k = —1) and to < 1/4, the exact solution for equation (IT51) is equal 
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where Cq is an integration constant. For the case of negative curvature (fc — —1) and m > 1/4, or 
positive curvature {k = +1) and m < 1/4, we can obtain the fohowing sohition for equation ([15]): 



a{t) = 



C 



cosh I 



where Cq is an arbitrary constant. It is interesting that in these cases, as it can be observed from 
the special feature of solutions and equation ([15]), this model behaves similarly to the FRW model 
in which the only source of gravity is the effective cosmological constant A = 3C/M. 

(II) To involve in our consideration the widely discussed modifications of Y-M theory (see, for 
example, [36]-[39], [l^), we assume that function ^(/) — B I™, and <&(/) arbitrarily depends on 
/. Then the main equations of our model can be written down as follows: 



(1 - 4m) 



2^ 

a 



-J <i 



4m) + — 



SttG [3$(/) -4/$'(/)] 



3B 



(16) 



a\2 k SvrG $(/) 

(1 - 4m) + 3 = ^ ' 



Combining these equations, it is possible to obtain, instead of ([1 
the second derivative of the scale factor: 

a _ SvrG [(1 + 2m)$(/) - 2 / $'(/)] 
3B 1 - 4m 



(17) 

the following equation for 
(18) 



It is easy to verify that (|T5)) is a differential consequence of equation ([rf]) which remains the only 
independent equation of our model. Nevertheless, from equation (IT^ it is possible to obtain a 
necessary condition for the accelerated regime (see Table 2). 



Ila 


B>0 


m > 1/4 


e(/) > 


lib 


B>0 


m < 1/4 


e(/) < 


lie 


B <0 


m > 1/4 


e(/) < 


lld 


B <0 


m < 1/4 


e(/) > 



Table 2. Conditions for m and function 8(/) — 2 /<&'(/) — (l + 2m)<i>(/), corresponding 
to accelerated expansion according to equation ([TS]) . 



Let us then consider the case of non-Abelian Lagrangian of the Born-Infeld type (BI) : 



f32 pa^paik [F-^F-^^Y 

^^^^=4^lV^^ 

where (3 is the critical intensity of Bl-field, F"^. is a dual Y-M tensor. From formulas ([4]) and metrics 
([2]), it follows that the second invariant of Y-M field for our solution Fi^F*'^ = 0. So, we can write 
down <&(/) as 

1 



where a = 1/2(3'^. From the latter, it is easy to find that 



,/ X / X 1 [(2m, + + 2al - 1) - 4m,a/l , , 

167ra ^1 + 2al 



In view of inequalities 2al > and y/1 + 2al > 1, one can find that expression (|20]) will be positive, 
only if m G (0, 1/4). Thus, as it follows from the explicit form of invariant ([7]), the positivity of 
expression (j^O)) will occur only after crossing the critical value of scale factor: 

"^')>"-=b(l-4m2) 
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One can see that this case corresponds to (lid) in Table 2, if B < 0. 

When the value of ([^0]) is negative, m remains arbitrary. However, the acceleration conditions 
are different for < 1/4 or > 1/4. The first case corresponds to (lib) with m g (—1/4, 1/4) 
and B > 0. The cosmic acceleration is possible while < a{t) < Ucr- In the case (lie) (i.e. for 
TO > 1/4, S < 0), the accelerated expansion is always possible. 

So non-trivial dependence of the model behavior on free parameter m requires a thorough 
research. Therefore, we are going to give more details and consequences of the model considered 
here in our further investigation. 

4 Conclusion 

In summary, the modified non-minimal Y-M theory in FRW non-flat cosmology are studied in 
this paper. First of all, we have derived the set of main equations which determines the model 
dynamics: ([6]), (l8|), ([9|). Throughout the last section of the paper, non-minimal coupling to gravity 
is described by the factor $(/) in the Einstcin-Hilbert sector of action ([1]). The non-trivial solution 
of the modified Y-M equation ([8|) proposed by one of the authors (V.K.S) earlier allows us to build 
several modifications of accelerated cosmic expansion in the frame of non-minimal coupling. Besides 
general study, we have considered in detail the power-law dependence of $(/) on its argument. We 
have derived the basic equations for the cosmic scale factor in our model, and have provided several 
examples of their solutions. This work implies that the cosmological applications of modified Y-M 
theory with non-minimal coupling to gravity may have more fruitful phenomena, which is worth 
studying further. 
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